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Abstract 

We derive the scattering amplitude N(r) for a QCD dipole on a dense target in the semi-classical 
approximation. We include the first subleading correction in the target thickness arising from ~ p 4 



operators in the effective action for the large- x valence charges. Our result for N(r) can be matched 



to a phenomenological proton fit by Albacete et al. over a broad range of dipole sizes r and provides 
a definite prediction for the ^-dependence for heavy-ion targets. We find a suppression of N(r) 
for finite A for dipole sizes a few times smaller than the inverse saturation scale, corresponding to 
a suppression of the classical bremsstrahlung tail. 
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In this paper we derive the dipole scattering amplitude N(r) [l| on a dense target in the 
semi-classical approximation. We restrict to a perturbative expansion of the Wilson lines 
valid at short distances r but include the first subleading (in density) correction arising from 
~ p 4 operators in the effective action for the large- a; valence charges. 

Our result may prove useful for a better theoretical understanding of the A-dependence 
of the initial condition for high-energy evolution of N(r), as discussed in more detail below. 

Let us first recall the setup for the McLerran-Venugopalan (MV) model {2]. The "valence" 
charges with large light-cone momentum fraction x are described as recoilless sources on the 
light cone with p a (x±,x~) the classical color charge density per unit transverse area and 
longitudinal length x~ . Kinetic terms for p are neglected since transverse momenta are 
assumed to be small. In the limit of a very high density of charge p the fluctuations of the 
color charge, by the central limit theorem, are described by a Gaussian effective action 

S MV [p]=jd^J_Jx . (1) 

Here, p?(x~)dx~ is the density of color sources per unit transverse area in the longitudinal 
slice between x~ and x~ + dx~ and 

K> 

dx'p\x-) ~ g 2 A x l* (2) 

X) 

is proportional to the thickness ~ A 1 / 3 of the target nucleus. This action can be used to 
evaluate the expectation value of the dipole operator 

D(r) = i-(tr V(x ± )V^y ± )) (3) 

9 A C F f j — 2 /„, — \ _2i__ 1 \ 



= ex P g^- J dx » ( x ) r lo S ) ^ 

= exp(-^!log^) , (5) 

where r = |xj_ — yj_ | and A is an infrared cutoff on the order of the inverse nucleon radius; the 
explicit r-dependence has been obtained in the limit logl/(rA) ^> 1. The scale Q s denotes 
the saturation momentum at the rapidity of the sources. For details of the calculation we 
refer to refs. and to the appendix below. 

To go beyond the limit of infinite valence charge density one considers a "random walk" 
in the space of SU(3) representations constructed from the direct product of a large number 
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of fundamental charges 6j. The effective action describing color fluctuations then involves 
a sum of the quadratic, cubic, and quartic Casimirs which can be written in the form [?J 

_ f p a (v±, v ± )p a (vi, v ± ) d abc p a (v^ } v ± )p b (v^, v_ L )p c (w 1 ", v ± ) 



s[ P ] = J d 2 v ± /°°^r{ : 



2fi 2 (v 1 ) K 3 

I 00 d P a (vi,v±)p a (vi, v ± )p b (vj, v±)/K, v±) , (0) 

-oo K i 



The coefficients of the higher dimensional operators are 

% ~ 9 3 A^ , (7) 
k 4 ~ <7 4 A , (8) 

and so involve higher powers of gA 1 ^ 3 . In what follows we restrict to leading order in l/«4 and 
drop the "odderon" operator ~ p 3 from (jSJ) since it does not contribute to the expectation 
value of the dipole operator at leading order in I//C3. The details of the calculation are 
shown in the appendix, here we just quote the final result for the T-matrix 



N(r) = 1-D(r) 



l2„2 1 r<2 -8 r poo 1 2 



log 



dz ji (z ) 



r 2 log 3 A- , (r 2 Q^<l). (9) 



4 rA Qir 3 k 4 

This now involves a new moment of the valence color charge distribution, namely 
J dx~ p, A (x~). We have calculated the 0(1/^4) correction analytically only in the "short 
distance" regime up to order ~ r 2 . Recall, however, that the effective theory (CQ) or ([6]) 
does not apply to the DGLAP regime at asymptotically short distances. Our result could 
in principle be extended into the saturation region r > 1/Q S by generating the color charge 
configurations p a (x~,xjj non-perturbatively, numerically [8]. 

The dipole scattering amplitude for a proton target has been fitted in ref. {9]] to deep- 
inelastic scattering data. The Albacete-Armesto-Milhano-Quiroga-Salgado (AAM QS) model 
for the initial condition for small-x evolution is given by 



N A AUQs(r,x = 0.01) = 1 - exp 



-I( r ^( IO ))^log(e + -y 



(10) 



with 7 ~ 1.119 . This model simultaneously provides a good description of charged hadron 



transverse momentum distributions in p + p collisions at 7 TeV center of mass energy 
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1 There are actually several fits, we refer to the AAMQS paper [9] for a more detailed discussion. 



This is a rather non-trivial cross-check: the MV model initial condition (J5]) overshoots the 
LHC data by roughly an order of magnitude at p± > 6 GeV lo| . 

However, since the model fllUp was introduced essentially "by hand" it is an open ques- 
tion how it extends to nuclei. This is a crucial issue for predicting the nuclear modification 
factor R p a ioi p + Pb collisions at LHC, and for heavy-ion structure functions which could be 



measured at a future electron-ion elC collider 



111 ] . One possibility is that the AAMQS mod- 



ification of the MV model dipole is due to some unknown ^-independent non-perturbative 
effect. Here, we explore another option, namely that for protons the effects due to the ~ p 4 
operators may not be negligible. 

We can match our result (|HD approximately to the AAMQS model by choosing 

1 2 ^ 



Cl g 8 

67T 3 Ql K 4 



dz /i 4 (z ) 



100 



(A 



(11) 



For nuclei, (3a ~ A~ 2 I 3 since each longitudinal integration over z~ is proportional to the 
thickness ~ A 1 / 3 of the nucleus while /i 2 (z~) and fx 4 (z~) are A- independent. The scattering 
amplitude for a dipole in the adjoint representation with this /3 is shown in fig. [IJ 2 One 
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FIG. 1: Left: scattering amplitude for an adjoint dipole (A^ = 2A — A 2 ) on a proton, 
assuming Q 2 = 0.168 GeV 2 and A 2 = 0.0576 GeV 2 . Right: same for a nucleus with 

A = 200 and Q 2 ~ A 1 ' 3 , A ~ A~ 2/3 . 



observes that the dipole scattering amplitude derived from the quartic action is similar 
to the AAMQS model over a broad range, rQ$ > 0.04. Discrepancies appear at very short 

2 To regularize the behavior at large r, for this figure we have replaced log -V — > log(e + i) and assumed 
exponentiation of the 0(r 2 ) expression. This does not affect the behavior at rQ s < 1. 
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distances where none of the above can be trusted. A more careful and quantitative matching 
of j3 to the AAMQS fit beyond the leading log 1/rA ^> 1 approximation should be performed 
in the future. 

On the right, we plot the scattering amplitude for a nucleus with A = 200 nucleons, 
assuming that Q 2 S ~ A 1 / 3 while 0a ~ A~ 2 / 3 . This illustrates how expectation values obtained 
with the quartic action converge to those from the MV model when the valence charge density 
is high (i.e., at large A 1 / 3 ). If our idea that the p 4 term in the action provides the explanation 
for the AAMQS model is indeed correct then their modification of the MV model should 
vanish like a ~ A~ 2 I 3 . This should be observable via R p a at the LHC. 
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I. APPENDIX 



Expectation values of operators 0[p] are computed as 

(0[p\) = Jv P 0[p] e~ s[p] I Jvp e~ s ^ . 
We work perturbatively in 1/^4 and keep only the first term in the expansion. Then, 



(0[ P ]) = 



JVp 0[p] e- s ^p\ 




JVp e- s Glp] 


1 


- £ / d 2 v ± f dv^dv^p a Vi p a Vi p b V2 p b V2 





(°[P] (l - £ / d ^± I dv^dv^p a Vi p a Vi p b V2 p b V2 ) } 
(1 " 5 I d?v L J dv^dv^p a Vi p a Vi p b V2 p b V2 ) G 
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In lattice regularization the denominator evaluates to 
1 



1 



1 - 



d z v ± I dv 1 dv 2 p a Vl p a Vl p b V2 p b V2 

K 4 J J I Q 



1 



(N 2 - If 



dv p (v ) 



+ 2(JV?-1)/ dv-p\v~) 



(13) 



«4 Avj_ 

where N s denotes the number of lattice sites (the volume) and Axj_ the transverse area of 
a single lattice site (square of lattice spacing). Also, we have used a local (pp) correlation 
function: 

(p a (x- } x ± )p b (y-,y ± )}=5 ab p 2 (x-)5(x--y-)5(x ± -y ± ) . (14) 



A. Dipole Operator 



We are interested in the expectation value of the dipole operator defined as 



D(x ± ,3/j.) = — tr y(x ± )0(y ± ) 



Here, V denotes a Wilson line 



V(x ± ) = Vexp^-ig 2 J dz-^p a (z-,x ± )t a ^ 



(15) 



(16) 



where 



If 1 

— r p a (z",xj.) = / c/ 2 z_ L G (x_ L - z±)p a (z~,z±) = — A + , (17) 

is proportional to the gauge potential in covariant gauge. The matrices t a are the generators 



of the fundamental representation of SU(3), normalized according to tr t a t b = \5 ab . 
Gq is the static propagator which inverts the 2-dimensional Laplacian: 

d 2 



2 G (x_l - z ± ) = 5(x ± - z ± ) ; 

d 2 k gik x .(x x -z ± ) 



dz^ 

G (x_l - z ± ) = - 



(2tt)2 k 2 ± 

With this propagator we can write the Wilson line as 



V(x_l) = Vexp!^-ig 2 J dz~ J d 2 z ± G (x ± - z ± ) p a (z~, z ± ) t a | 



(18) 
(19) 

(20) 



The correlator (V(xj_)V"''(yi)) for a Gaussian (MV) action has already been calculated 
before, see for example ref. la]. The result is: 



(V(x ± )V\y ± )) G = e W \-^(tH a ) 



{-(a) dz-p 2 (z~) y dV [G (x ± - z ± ) - G (yx - z ± )] 2 | 

(21) 



Note that this is diagonal in color (proportional to 13x3)- To calculate the expectation value 
of the dipole operator with the new action, we first expand the Wilson lines order by order 
in the gauge coupling g 2 , 

V(x ± ) = 1-tg 2 [ d 2 z ±1 G (x ± - z±1 ) / dzip a ( Zl )t a 

J J — CO 

/pCO nOO 
d 2 z ±1 d 2 z ±2 G ( X± - z ±1 )G (x ± - z ±2 ) / dz-{ / dz 2 p a ( Zl )p b (z 2 )t a t b 
J— 00 J z 1 

+ ■■■ (22) 
\/t( y± ) = 1 + ig 2 f d 2 u L1 G (y ± -u L1 ) duip a ( Ul )t a 

J J —OO 

-# 4 [ d 2 u ±1 d 2 u ±2 G (y± -u ±1 )G (y± - u ±2 ) / duf / du 2 p a (u l )p b (u 2 )t a t b 
J J —00 J —00 

+ ■■■ (23) 

For brevity we only write the terms up to 0(g 4 ) but below we shall actually require terms 
up to 0(g 8 ). 

To zeroth order in g, the expectation value of the correlator is just 1: 

O(g )=l . 
The order g 2 contribution is zero because (p a (z)) = 0: 

O(g 2 )=0 . 

B. Order g A 

The first non-trivial contribution arises at 0(g 4 ) and is given by the sum of expectation 
values of three terms (two from each Wilson line and one mixed term): 

/poo nOO 
d 2 z ±1 d 2 z ±2 G {x± - z ± i)G (x± - z ±2 ) / dz^ dz^ (p a ( Zl )p b (z 2 ))t a t b , (24) 
J — 00 ./ z x 

-g A d 2 u ±1 d 2 u ±2 G (y ± -u ±1 )G (y ± -u ±2 ) / du^ / d«J (p a ( Ul )p b (u 2 ))t a t b , (25) 

J J —CO J —CO 

/n pco pco 

d 2 z ±1 / d 2 u ±1 G (x ± -z ±1 )G (y ± -u ±1 ) / dz^ / du^ {p a ( Zl )p b ( Ul ))t a t b . (26) 
J J —CO J —CO 
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Using (|T2|) and (f2~T]) the first term becomes (we will divide by the normalization factor 
later) : 

- dz-fi 2 (z-) J rfV^(xi-zi) 

4 (■ 

+ — d 2 z ±1 d 2 z ±2 d 2 v±G (ic± - z ± i)Go(x_l - z ±2 ) 

«4 J 



tfar / rf^ / dv^dv^{ptoloto c v o d v o d Jt a t b 



All possible contractions for the second term in the above expression are 

(ptAAAAAJ = (pIA) HMipiM + 2 (ptMiPnPi)} 
+4 ( P a z M [( P b Z2 P c J(pip d J + 2 (p b yj(p c vl p d J] , 



(27) 



(28) 



and are shown diagrammatically in fig. [21 Using (fl4l) to perform the contractions and 



o u 




0> 1,2 ^ O" 

FIG. 2: (V) at order c/ 4 //€ 4 . 



dividing also by the normalization factor (|T3l) leads us to 3 

1 



i jv s 



«4 



£t ( N ? - !) 2 Too dv-^(v-) + 2 (JV* - 1) dv-^{v~) 



1 - 



1 iV s 



^ 4 t°ta 



^4 Av_i_ 
J d 2 z ± G 2 {x ± - z x ) 



(N 2 - If 



00 n 2 

dv~ jji 2 (v~) 



+ 2 (iV c 2 - 1) / dv-fi\v~) 



(N 2 -l) / dz~\i 2 {z~) / *;>>-) + 2 / dz-fi 6 {z-) 



00 
00 



K4 Avi 

The third line in the above expression cancels 4 the normalization factor once the latter is 



+ 



(29) 



3 One has to be careful with performing the integration over z^'. dz^b{z\ — z%) = 1/2. If x~ is 
discretized, should be placed ahead of z-f by at least half a lattice spacing Ax~ . Similarly, when 
expanding a single Wilson lines to order g 6 : z^ > z^ + Ax~ /2 > z-f + Aa; - . 

4 This is the standard cancellation of disconnected diagrams. 
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expanded to leading order in I//C4 so that the previous expression simplifies to 

- 1-*%/ ^(xi-zj 

dz~ p?(z~) 



4 



(iV c 2 -l) / dz~p 2 {z~) / ^"/i 4 (t;-) + 2 / dz-/i 6 (z-) 



.(30) 



The correction is absorbed into a renormalization of J™ dz p 2 (z ) in order that the two- 
point function (pp) remains unaffected by the quartic term in the action: 



dz p 2 (z ) 
dz~ p 2 (z~) — 



(N 2 -l) / dz-p 2 (z-) / dv-p\v-) + 2 / dz-p 6 (z-§ol) 



K 4 Avj_ 

Finally, the expectation value of V(x±) to order g A can simply be written as 



Similarly, 



-ytX J°° dz-~p 2 (z-) J dV^( Xl - Zi ) . 
- ytX dz-p 2 (z-) J d 2 u ± G 2 (y ± - u ± ) . 



(32) 



(33) 



The mixed term will be the same as the previous terms except with a positive sign and 
without the factor of 1/2 which originated from the path ordering (z-f and are not 
ordered relative to each other): 

(34) 



gH a t a / dz-j?(z-) / d 2 z ± Go(x ± -zjGo(y±-z ± ) . 
Summing ( 1321) . ( 1331) and ( 1341) . we obtain the complete result at order g 4 : 

0(g A ) = -^t a t a J°° dz-~p 2 (z-) J d 2 z ± [Go(x ± - z ± ) - G (y ± - z ± )f 



(35) 



This is identical to the result obtained in the Gaussian theory once the two-point function 
{pp) ~ p 2 has been matched. This was to be expected, of course, since only two-point 
functions of p arise at 0(g 4 ). Note, also, that fl35|) vanishes as y_i_ — > x_l. 
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FIG. 3: Order g 8 contribution to (V). 

C. Order g 8 

Next, we consider order g 8 . There are all in all five terms: two of order g 8 from the 
expansion of a single Wilson line, two mixed terms (order g 2 from the first line and g 6 from 
the second line and vice versa), and one term from multiplying g 4 terms from both Wilson 
lines. 



1. g 8 from ^(xj_) 



First, we calculate (V^xjJ) at order g 8 . In the Gaussian theory 



77 (t a t a ) 2 



dz p?{z ) 



J d 2 z±Gl(x± - z. 



Again, using ([12]) . the correction is 



X 



d 2 z ±1 d 2 z ±2 d 2 z_ L3 d 2 z_ L4 / d 2 v±G (x± - z ±1 )G Q (x± - z ±2 )Go(x± - z_l 3 )G (x ± - z ±4 ) 



dz x I dz 2 



dzn 



dz' A 



dv- x \ ^ 2 -(p: l p^= 3 pf 4 p: l p^ lP 4p4)n 6 t c ^(36) 



All possible contractions for (p^P^P^pt.P^P^P^P^) are shown diagrammatically in fig. |3j 
The disconnected diagrams [3a] give the following contribution: 



7T (^) 2 



1 N. 



dz p}{z ) 



J d 2 z ± Gl(x ± - 





" POO 


{Nl - If 


/ dv~p, 2 (v~) 




_J — oo 



K 4 Avj_ 

which will cancel the normalization factor. 



+ 2(iV c 2 -l)/ dv-p\v-) 



(37) 
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The connected diagrams from fig. l3bl again renormalize fi 2 (z ) to p, 2 { 



dz jj, (z ) 



(N 2 - 1) 



dz fJL (z ) 



dz jj (z ) 



Finally, the diagrams [3c] give the correction 

9 S 



dz~ fi 4 (z~) / rf 2 z_|_GQ(xj_ — z_i 
) J J 

That leads us to the complete 0(g 8 ) contribution to (V(x±)): 



- - (t a t a ) 
^4 



dz fl (z ) 



9_ 



00 
00 



dz fl (z ) 



J d 2 Z_LGo(x_L - Zj 



(40) 



The expectation value of V\y±) is obtained from the above by substituting xj_ — > y±. 



dv-jj^v-) + 2 / dz-fi 2 (z~) / dv~fi b (v-j-\8) 



(39) 



2. g 6 from V(x_\_) x g 2 from V^(y±) 

Next is the mixed term obtained when multiplying the g 6 term from the x Wilson line 
with the g 2 term from the y Wilson line. In the Gaussian theory, 



-y(tX) 2 



The correction is: 

,8 



-1 2 



dz fi 2 (z ) 



J d 2 z ± G 2 (x ± -z ± ) J rfVG (x ± -u ± )G (y ± -u ± ). (41) 



+ ~ J d 2 z ±1 d 2 z_ L2 d 2 z_ L3 G (xL± - z±i)G {x± - z_ L2 )G , (x_ L - z± 3 ) J d 2 u ±1 G (y_L - u ±1 ) 

poo poo poo poo 

dzi / dzj / dz 3 I dui I dvi 



X dV, 



dVn 



x (plAAXAA^LpOt-thH* . (42) 

The possible contractions are given in fig. HI As before, the disconnected diagrams Ha 
cancel against the normalization, while the diagrams in 14b I renormalize fi 2 to ft 2 ; finally the 
diagrams He] will give the correction. Summing all these diagrams plus the Gaussian part, 
we get: 



y it a t a f 



dz fl (z ) 
,8 



+4 9 - {t a t a f 
K4 



J d 2 z ± Gl(x±_ - z ± ) J rf 2 u_LGo(x± - u_L)G (y_L - u ± ) 

dz-fi\z-) I / d 2 z ± G 3 (x ± - z ± )G (y± - z ± ) (.43) 
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FIG. 4: V r (x_ L )y t (y_ L ) at order g 8 / ' k a (g e from V(x_l) x g 2 from V^(y±)). 

Once again, a similar contribution (with xj_ -H- yj_) arises from V(x_|_) at 0(g 2 ) times V^(y±) 
at order g 6 . 



3. g 4 from V(~x±) x g 4 from V'(y±) 

The last term to consider is the one obtained when multiplying 0(g 4 ) from the x Wilson 
line with 0(g 4 ) from the y Wilson line. The Gaussian contribution is: 



+y (t a t a y 



dz ji (z ) 



dz fi 2 (z ) 



J d 2 z±Gl(x± - z ± ) y rf 2 u ± Go(y ± - u ± 
J d 2 z ± G (x 1 _ - z ± )G (y_L - u ± 



(44) 



12 



D 2 «1 





(a) 



Vi Vi 
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FIG. 5: Expectation value of V(xj_) at order g A times V(y±) at order g 4 



For the quartic action we have to add 



9_ 



J d 2 z ±1 d 2 z ±2 G (x ± -z_ l1 )G (x_l -z_l 2 ) J d 2 u ±1 d 2 u ±2 G (y ± - uj_i)G (yj. - u ±2 ) 



X / C? 2 V_L / 



dw 1 



du 1 



dv 



X \P Zl Pz 2 Pui Puo Pvi Pvi Pvo Pvn / * * t 1 ■ 



(45) 



Following the same procedure as before, the diagrams from fig. [5] give: 



\ (t a Q 2 



-i 2 



dz fx (z ) 



-6 9 - (t a t a ) 
K4 



-1 2 



// 2 (;z ) 
c/z~/i 4 (z~) 



n 2 



rf 2 z±G (x± - z_L)G (y_L - uj_) 



/ d 2 z ± G 2 (x ± - z ± )G 2 ( yi _ - z ± ) . 



(46) 
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4- Complete order g' 



Combining eqs. (HOlfx i <-> yj_), (1431-1- x i -B- and (|46|) we get the complete contribution 
to the expectation value of the dipole at order g 8 : 



2 V 2 



-- {t a t a f 

K4 



dz~fi 4 (z~) 



J d 2 z ± [G (x ± - zl) - G (y± - z ± )] 2 
y rf 2 z ± [G (x ± - z ± ) - G (y ± - zj]' 



Note that this again vanishes as y_i_ — > Xj_. 



D. Complete expectation value of the dipole operator 

Adding together the terms of order 1, g 4 and g 8 , the expectation value of the dipole 
operator becomes 

D(x x -y ± ) = (^-tr V(x ± )0(y ± )) = 

1 - \c F dz-f{z-) J rf 2 z ± [G (x ± -z ± )-G (y ± -z ± )] 2 



1 (g 4 C F ^ 2 



00 1 2 r 

r.2 



dz jl (z ) 



t 

K4 



dz ji (z ) 



1 d 2 z ± [G (x ± - z ± ) - G (y ± - z ± )f 
J d 2 z ± [G (x ± - z ± ) - G (y± - z ± )] 4 + • • ■ (47) 



where 



iv 2 -i 

2iV„ 



To write this in a more compact form we introduce the saturation scale 

2tt' 

so that 



9 

Q 2 = —Cf I dz~jl 2 (z~' 



(48) 



D(r) = l-nQ 2 J d 2 z ± [G (x± - z ± ) - G (y± - z ± )f 
y rf 2 z ± [G (x ± - z ± ) - G (y± - z ± )] s 



a 8 

- g -c 2 F 

K4 

Here, r = |xj_ — yj_|. 



dz \i 4 {z 



j d 2 z ± [G (x ± - z ± ) - G (y ± - z ± )f 



(49) 
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E. Explicit evaluation of D(r) to leading logl/r accuracy 



It is useful to obtain an explicit expression for D(r) in the limit logl/r A ^> 1, where A 
is an infrared cutoff on the order of the inverse nucleon radius. 
The first non-trivial term in eq. f)49p gives 



ttQ 2 J d 2 z ± [G (x ± - z ± ) - G (y± - z ± )f 

oo 

«/*[l-W]-^log-L, 



(50) 



in the leading log 1/rA ^> 1 approximation. 
Next, we need to compute the integral 



J d 2 z ± [G (x ± - z ± ) - G (y± - z±)] 4 • 



From eq. (fT5|) for the propagator, 



J d 2 z_LGo(x_L - z ± ) = j^y, J d 2 z± J (fk^^d^cfki ^2^2^.2 6 



1 



(2tt)< 



d 2 Kd 2 kd 2 Qd 2 q 



6(K + Q) 



(27T) 



d 2 K 



(f + k)^(f-k)^(f + q )^(?- q ) 
1 



(51) 

t(ki +k 2 +k 3 +k 4 ) ■ (x_l 

(53) 



-1 2 



d z k- 



(f + k) 2 (f-k) 



(54) 



We regularize the integral in the square brackets by introducing a cutoff A: 

1 2tt 



d z k- 



Then, 



(f + k) 2 + A 2 (f-k) 2 + A 2 



d 2 K 



d 2 z ± G Cx ± - z ± ) 



K 2 ^ 



K 2 



K 2 

A 2 



(55) 



(2tt) 4 J K 4 A 2 
dK . o K 2 



1 1 



(2tt) 3 7 a iT 3 A 2 (2vr) 3 A 2 



(56) 



Following the same procedure for J d 2 z±Gl('Xj_ — z±)G 2 ] (y± — z±) we arrive at: 

/ d 2 z ± G 2 (x ± - z ± )G 2 (y ± - z± ) = JL I ^e^-(-) log 2 



K 2 

A 2 " 



(2vr) 3 J A K 3 



00 c? 2 K 2 K 2 1/1 1 o „ s 

Jo(Kr) log 2 — = j— rr + -r 2 log 3 (rA) 



A 2 (2tt) 3 VA 2 3 



[r 2 log 2 (rA)j 
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Similarly, 



J d 2 z ± Gl(x ± - z ± )G (y± - z ± ) 



(27T) 



d 2 K 



K 2 



if -(xx-yx) j 

K 2 g A 2 



e «q-(xx-y±) 



(f + q) 2 (f-q) 



Using: 



ci 2 q 



giq-(xx-yx) 



(I + q) (f-q) 

1 + iq ■ r 



we get 



/ 



<Pz ± Gl(x± - Zx)G (y± - z ± ) 



1 



(2tt) 3 
so that, finally, 

D(r) 



oo t2 



rf 2 K 



K 2 



Jn -At = - 



1 1 



A 2 4(2vr) 3 VA 2 3 



4 1 

— + -r 2 log 3 (rA) 



r 2 2 
1-^log— + 



1 1 9 8 „2 



rA 6tt 3 K4 



Ci 



dz fi A (z ) 



r 2 log 3 



rA 



1 _ i g _L + /3 r 2g2 log s i_ 



(57) 



(58) 



(59) 

(60) 
(61) 



/S has been defined in eq. (jTTj) . 

Performing a Fourier transform we obtain the transverse momentum dependence of the 
dipole for fci > Q s > A: 



D(k A 



.01 



2ir ^± - J—C 2 



dz fi 4 (z ) 



1 h 2 

fcT log A 2 



(62) 
(63) 



The first term was taken from appendix B of ref. The second term provides a correction 
to the classical bremsstrahlung tail for finite valence parton density. 



F. Gluon density 



One can define 
small Xj_, as 



4j a Weizsacker- Williams like gluon density of a nucleus, in the limit of 

(64) 
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Z2 Z 3 




FIG. 6: (V) at order g / K4 for the action ([6 
For the quartic action 

N(xl) = ^ log -L - /^Q 2 log 3 -L . (65) 
The first term, via eq. flSljl . gives xGyi ~ A [4|. The second term gives a correction 

G. Form of the quartic term in the action 

In this section we explain why the p's in the quartic term of the action should sit at two 
different points in the longitudinal direction in order that N(r) vanishes as ~ r 2 as required 
by color transparency. 

First, let us note that the correction to the dipole scattering amplitude due to the quartic 
term in the action results from the diagrams |3c] He] and [5c] Now, let us consider the action: 

P a (v~, vi_)p a (v~, V_l) 



S[p] = J d 2 v±J dv- |- 



l P a (v ,v±)p a (v ,v±)p b (v ,v±)p b (v ,v±) | ^ 

The analogue of diagram [3c] for (V) at order g 8 for this action is shown in fig. [6] However, 
this diagram vanishes due to the longitudinal path ordering in the Wilson line. 

The same reasoning applies to the analogue of He] shown in fig. [7] In this case the three 
points Zi , Z2 and z% , can not be connected simultaniously to one v~ point. The delta 
functions coming from this kind of contractions, 5(zi — v~)5(z2 — v~)5(z^ — v~), imply 
overlap of z± and z% which is not there in a path ordered integral. 

That leaves us with only one type of diagram proportional to I//C4, shown in fig. [8] This 
diagram is not zero since there is no relative ordering between the points z~ and u~ on the 
two lines. This diagram has a constant r-independent contribution which does not cancel 
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"1 




FIG. 7: Order g 6 from V^xjJ times order g 2 from V'(y_i_) at order I/K4. 




FIG. 8: Order g 4 from V(xj_) times order g 4 from V^yjJ at order 1/k 4 . 

because of the missing diagrams |6] and [3 This is (only terms at order l//c 4 are given) 
D{r) „ l-^l(C 2 F - A ) r dz-p\z-) [ rfV^(x ± - Z JG 2 (y ± - Z± ) 

Hence we see that for the action fl67|) that N(r) = 1 —D(r) approaches a constant as r — > 0, 
in violation of color transparency. 

For numerical (lattice gauge) computations of expectation values ( |T2i) it may be easier 
to integrate p{x~) and to drop the longitudinal path ordering, 

00 

«*x) S /<fa-/<*-,*J (69) 

— OO 

\Z(x_l) = e v i . (70) 



For a detailed discussion, we refer to ref. 
action 



. One would then consider the two-dimensional 



M= f^I^^ + ^E^iM£M\ . (71) 

J [ 2/i K4 J 

The diagrams from figs. [6] and [7] then do exist and cancel the r-independent contribution 
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from fig. [8] so that again N(r) ~ r 2 at r — > 0. 
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